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Wilson loops in large N gauge theory exhibit a weak to strong coupling transition as the loop is dilated. A multiplicative 
matrix model captures the universal behavior associated with this transition. A universal scaling function is obtained 
in a double scaling limit. Numerical studies show that both large N QCD in three dimensions and the SU(N) principal 
chiral model in two dimensions are in the same universality class. 
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Large N gauge theories may provide a bridge between QCD and string theory. We hope that this would produce 
a way to make a calculational connection between some weak and strong coupling phenomena of continuum QCD. 
Recent progress [H, E|, [^l in this direction is the main subject of this talk. 

We will consider large N gauge theories on a d = 2,3,4 dimensional continuum torus. It is sufficient for our 
discussion to consider a symmetric torus of length I in all four directions. We will consider the 't Hooft limit [4[ 
with a finite number of fundamental fermion flavors. There is no back reaction from the fermions on the gauge field 
dynamics. 

In all three dimensions, there is a critical size l\ such that continuum reduction [5[ holds for I > l\ and there are no 

finite volume effects. The critical size, l\, is equal to T c where T c is the deconfining temperature. Two dimensions 

is a special case since l\ = and large N QCD is always in the confined phase and there is no dependence on the 

temperature. Therefore, continuum Eguchi-Kawai reduction 0] holds in two dimensions 0. We will focus on the 

confined phase of large N QCD in this talk. 

Let us first make a remark about the Gross- Witten transition [8] in two dimensional large N lattice gauge theory. 

This transition has analogs in large N lattice gauge theories in three and four dimensions. The eigenvalue distribution 

of the plaquette operator opens up a gap around —1 at the transition point enabling the presence of topologically 

disconnected classes of gauge fields in the continuum limit Q . These transitions do not have a continuum analog. 

We will assume that we are on the physical side of this transition when discussing large N lattice gauge theories. 
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The stage is now set to discuss the transition from weak coupling to strong coupling in the confined phase of 
large N QCD. We will use the Wilson loop operator, W, a unitary operator for SU(N) gauge theories, to probe the 
transition. Large (area) Wilson loops are non-perturbative and correspond to strong coupling. Small (area) Wilson 
loops are perturbative and correspond to weak coupling. 



(N 
> 



X 

5-H 



The probe is defined as 

W N (z,b,L) = (det(z-W)) (1) 

and is the characteristic polynomial associated with the operator. W is the Wilson loop operator, z is a complex 
number, N is the number of colors, b — is the lattice gauge coupling and L is the linear size of the square loop. 
(■ ■ •) is the average over all gauge fields with the standard gauge action. 

The eigenvalues of W are gauge invariant and so is the characteristic polynomial. The eigenvalues lie on the unit 
circle and all of them will be close to unity for small loops. The eigenvalues will spread uniformly over the unit circle 
for large loops. The characteristic polynomial exhibits a transition at N = oo when L — > L c (b). This is a physical 
transition since L c (b) will scale properly with the coupling, b, as one approaches the continuum limit. 
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The eigenvalue distribution of the Wilson loop operator in two dimensional QCD was shown to undergo the above 
transition by Durhuus and Olesen and further relevant discussions can be found in 11 , ljj • The scaling function 
in the double scaling limit was derived for two dimensional large N QCD in Q. Two dimensional gauge theory 
on an infinite lattice can be gauge fixed so that the only variables are the individual plaquettes and these will be 
independently and identically distributed. The Wilson loop operator, W, can be written as W = Ylj=i Uj where the 
Uj's are the transporters around the individual plaquettes that make up the loop and n = L 2 is equal to the area 
of the loop. The measure associated with any Uj can be set to P{Uj) = J\fe~~ Tr H j where Uj = e leHj and e plays 
the role of gauge coupling. The dimensionless area is given b y t = e 2 n which is kept fixed as one takes the limit 
n — ► oo and e — > 0. This is the multiplicative matrix model of j 1 31 ] - In the continuum limit, the parameters b and L 
get replaced by one parameter, which is denoted by t in the model, and the characteristic polynomial becomes 



W N (z,b, L)^Q N {z,t) 
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A rearrangement of Q^(z,t) for SU(N), 



Z N (z, t) = Q N (z, t)(-l) N e^ (-*)"* = elni ~ Z) E ' El>J ' W 

CT1,(T2,...<TJV = ±5 

puts it in the form of the partition function for a spin model with a ferromagnetic interaction for positive r with 



ln{— z) as a complex external magnetic field. Therefore, the conditions for Lee- Yang theorem 1J] are fulfilled and 
all roots of Qn( z , t) lie on the unit circle for SU(N). 

The transition from weak coupling to strong coupling can be intuitively seen using the characteristic polynomial, 
Qisr(z,t). In the weak coupling (small area) limit we have t — and Qx(z,t) = (z — 1) N . Therefore, all roots are 
at z = 1 on the unit circle. In the strong coupling (large area) limit we have t = oo and Qn(z,£) — z N + (— 1) N . 
Therefore, all roots are uniformly distributed on the unit circle. 

Qw(z,t) is analytic in z for all t at finite N. But, this is not the case as N — * oo and leads to a transition from 
weak to strong coupling in the N — * oo limit. 

There is a critical area, t = 4, such that the distribution of zeros of Qoo(z,t) on the unit circle has a gap around 



z = — 1 for t < 4 and has no gap for t > 4 lCJ, [l3|. To see this, we note that the integral representation ([2]) is 



dominated by the saddle point, v = A(t, z), given by 

A = A(M)= * - (5) 

With z = e 10 and w = 2A + 1, p(9) — — ^Re w gives the distribution of the eigenvalues of W on the unit circle. 
The saddle point equation at z = — 1 is 

w = tanh —w (6) 
4 

showing that w admits non-zero real solutions for t > 4. 

As N — » oo, one can define a scaling region around t = 4 and z = — 1 by 

*=TtV; z = ~e^ (7) 
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a and £ are the scaling variables that blow up the region near t = 4 and z = — 1. We can show that 

/4iV\* AT (JV-l)r N f°° 4 2,, 

Im^i^—j (-l) N e L ^(~z)-^Q N (z,t)= j due~ u ~ au +*» = C(£, a) (8) 

which is the scaling function in the double scaling limit associated with the characteristic polynomial. 

This scaling behavior is expected to be universal and was numerically tested using the lattice formulation in 
three dimensional large N QCD 0] and in the two dimensional SU(N) principal chiral model The precise 
statement of the continuum large N universality hypothesis is as follows. Let C be a closed non-intersecting loop: 
Xfj,(s),s G [0,1]. Let C(m) be a whole family of loops obtained by dilation: x^(s,m) = — x^(s),with m > 0. Let 
W(m,C(*)) = W(C(m)) be the family of operators associated with the family of loops denoted by C(*) where m 
labels one member in the family. Define 

N (y,m,C(*)) = (det(e= + e _ *W(m,C(*))) (9) 

Then our hypothesis is 

Even though the hypothesis has not been explicitly tested in four dimensional large N QCD, there is enough 
numerical evidence in suppport of the hypothesis Assuming the universal behavior, one can address the question 
of matching the universal data in the transition region to the perturbative side. 

Acknowledgments 

R.N. acknowledge partial support by the NSF under grant number PHY-055375 at Florida International University. 
H. N. acknowledges partial support by the DOE, grant # DE-FG02-01ER41165, and the SAS of Rutgers University. 

References 

[1] R. Narayanan and H. Neuberger, JHEP 0603, 064 (2006) |arXiv:hep-th/0601210| . 
[2] R. Narayanan and H. Neuberger, JHEP 0712, 066 (2007) |arXiv:0711.455ll [hep-th]]. 
[3] R. Narayanan, H. Neuberger and E. Vicari, JHEP 0804, 094 (2008) [arXiv:0803. 3833 [hep-th]]. 
[4] G. 't Hooft, "A PLANAR DIAGRAM THEORY FOR STRONG INTERACTIONS," Nucl. Phys. B 72, 461 
(1974). 

[5] R. Narayanan and H. Neuberger, Phys. Rev. Lett. 91, 081601 (2003) [arXiv:hep-lat /03030231 . 
[6] T. Eguchi and H. Kawai, "Reduction Of Dynamical Degrees Of Freedom In The Large N Gauge Theory," Phys. 
Rev. Lett. 48, 1063 (1982). 

[7] G. Bhanot, U. M. Heller and H. Neuberger, "The Quenched Eguchi-Kawai Model," Phys. Lett. B 113, 47 (1982). 
[8] D. J. Gross and E. Witten, "Possible Third Order Phase Transition In The Large N Lattice Gauge Theory," 

Phys. Rev. D 21, 446 (1980). 
[9] J. Kiskis, R. Narayanan and H. Neuberger, Phys. Rev. D 66, 025019 (2002) [arXiv:hep-lat /0203005] . 
[10] B. Durhuus and P. Olesen, Nucl. Phys. B 184, 461 (1981). 

[11] M. R. Douglas and V. A. Kazakov, Phys. Lett. B 319, 219 (1993) [arXiv:hep-th/9305047] . 

[12] A. Bassetto, L. Griguolo and F. Vian, "Instanton contributions to Wilson loops with general winding number 

in two dimensions and the spectral density," Nucl. Phys. B 559, 563 (1999) |arXiv:hep-th/9906125|. 
[13] R. A. Janik, W. Wieczorek, J. Phys. A: Math. Gen. 37, 6521 (2004). 
[14] T. D. Lee and C. N. Yang, Phys. Rev. 87, 410 (1952). 



'3Na 2 {C(*)) 



C(£,<*) 



(10) 



3 



